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Abstract 

Multiparticle boson states, proposed recently for 'independently' 
emitted pions in heavy ion collisions, are reconsidered in standard 
second quantized formalism and shown to emerge from a simplistic 
chaotic current dynamics. Compact equations relate the density oper- 
ator, the generating functional of multiparticle counts, and the correla- 
tor of the external current to each other. 'Bose-Einstein-condensation' 
is related to the external pulse. A quantum master equation is advo- 
cated for future Monte-Carlo simulations. 







1 Introduction 



Multiparticle production in high energy particle collisions is dominated by 
classical statistics. Bose-Einstein statistics of pions, nevertheless, proved to 
lead to quantum coherence effects which survive in the final multiparticle 
states [0]. The corresponding mechanism can be handled within standard 
quantum statistics. Yet, some recent works |fj, |3| continued, mainly for 
historic reasons as Weiner || points out, to use the "traditional" method 
of wave function with tedious explicit Bose-Einstein-symmetrization. We 
shall see how standard second quantization methods lead to the correct re- 
sults in a shorter and clearer way. I shall discuss the dynamical conditions 
of Bose-Einsein condensation, and I outline a master equation suitable for 
Monte-Carlo simulations. 

For the concentrated study of effects of Bose-statistics, a simple scheme of 
'independent' emission has been proposed f|, || H, [7j. Sec. 2 recapitulates the 
features of these multiparticle states, sparing the burden of separate 'sym- 
metrization'. Sec. 3 introduces the generating functional of the measured 
counts H directly in second quantized formalism. Introducing chaotic clas- 
sical currents, advocated e.g. by Ref.[|l[] and used e.g. in [^, |10[ , I construct 
the simplest quantum dynamics reproducing the corresponding multiparti- 
cle quantum states (Sec. 4). The measured multidetector counts turn out 
to be identical to the corresponding spectral intensities of the effective cur- 
rent (Sec. 5). I show that the existence of Bose-Einstein condensate imposes 
explicit analytic constraints on the intensity and on the spectrum of the 
external effective current (Sec. 6). Finally, I generalize the simple quantum 
dynamics and propose a quantum master equation suitable to efficient Monte 
Carlo simulation of the multiparticle density matrix itself (Sec. 7). The Letter 
concludes with summary (Sec.8). 

2 Independent multiboson states 

When searching for a class of multiparticle density operators p represent- 
ing independent bosons, consider first the Gibbs canonical state for non- 
interacting bosons at inverse temperature /3. The bosons remain independent 
if, formally, we assign different instantaneous temperatures /3k to each modes 
k, i.e. we assume p ~ ( — yj k A^ak) . Moreover, the bosons remain indepen- 
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dent if we assign different temperatures to a generic (maybe non-stationary) 
set of orthogonal modes instead of the momentum eigenstates. Hence, we 
arrive at the following class of 'independent multiboson states' (IMS): 



p = det(l - e /3 ) exp - ^ Ac'^'Ok (1) 

V k,k' / 

where f3 positive matrix. Let us define the 1-particle density matrix from the 
above state: 

(g) 

Pk'k - „ 7j | v • {*) 

Ek(k|p|k) 

Using Eq. (1), we find the following matrix relation: 

e'l 3 = vp {v = tre~ p ) . (3) 

The IMS (1) can be rewritten in terms of the 1-particle density matrix p and 
the parameter v (whose physical interpretation remains a bit involved): 



p = det(l - e up ) exp ( v ^ pk^a^ <g> a k ] p • (4) 

k,k' 



This form might give an insight into the kinematics of the particle creation 
from the vacuum p . 

We have to note that the IMS are non- stationary quantum states. Yet, 
the measured quantities = a\xik are not sensitive to the time evolution of 
the IMS p. This will be formulated in the next Section. 



3 Generating functional 

We introduce generating functionals for the multiparticle final state momen- 
tum distributions. A compact heuristic form of definition is the following: 

G[u\ = tr (pnW k ) , (5) 

where are auxiliary variables and = a^a^. If we introduce the diagonal 
matrix u by Wk'k = ^k'k^k then, using the IMS density operator (4), the 
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generating functional takes the following form: 



G[u] = P£ ~ VP \ . (6) 
L J det(l - vup) v ; 

The logarithmic generating functional g = logG can be expressed through 
its Taylor-expansion in a transparent way |5], [?|: 



9M = E - (tr(up) r ~ trp r ) . (7) 

r=l r 

For «k = u, it yields the (logarithmic) multiplicity generating function 

oo r 

9(u) = E -trp r (u r - 1) (8) 

whose Taylor-coefficients are the combinants (c.f. [|TT|j ). 

The derivatives of the generating functionals at u = yield the exclusive 
distribution/correlation functions. In experiments, we can easily measure 
the inclusive distributions instead, which are the derivatives at u — 1 f|. To 
make these derivations more convenient, let us substitute vp in the generating 
functionals (6-8) by vp = a/(l + a), where a will be the correlation matrix 
of 'currents' of Sec. 4: 

G\u] = — — i r --t , (9) 

1 J det[l -{u- l)a)\ V ; 

oo v r 

g[u] = — tr(ua — a) r , (10) 
T^i r 

oo r 

g(u) = j:-trp r (u-iy . (11) 

r=l r 

Comparing these expressions with the Eqs. (6-8) we see that the inclusive 
distributions/correlations will depend on the current correlation matrix a 
exactly the same way as the exclusive distributions/correlations depend on 
(z/-times) the 1-particle density matrix p. 
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4 Dynamics 



Multiparticle production can not be derived from 'first principles'. 1 can 
certainly not overcome the well-known difficulties. Instead, I present the 
simplest quantum dynamics which produces exactly the class (I) of IMS. I 
postulate the following effective Hamiltonian: 

H = 5(t)J2(-J^ k + J k a£) , (12) 

k 

where J k denote the Fourier-components of a certain effective external field 
J exciting the boson-field. The 'current' J(x) = J(x)<5(£) is non-zero in the 
collision area and we assume that the collision time can be taken infinite 
short. Let us calculate the unitary effect of the Hamiltonian (5) on the 
vacuum: 

Wj) = exp J k«k + Jk4)) |0) = exp (^El J k| 2 -*E J k4) |0> , 

(13) 

which is otherwise a product coherent state Ilk® I — zJ k ). 

These final states are pure states whereas the IMS are mixed ones. 
Obviously, no unitary dynamics can create mixed states from pure ones. 
Therefore, I consider unitary dynamics in random external fields: I assume 
Gaussian distribution for the stochastic fluctuations of the current J. Let the 
mean values M[J k ] be always zero. Also we assume that M[J k / J k ] = 0, which 
is equivalent to a random phase for all J k . We denote the only non-vanishing 
correlations by the non- negative Hermitian matrix a: 

M[J k ,JZ] = ov k . (14) 

After these preparations, we can define the density operator p of the final 
state as the stochastic mean value of the pure coherent states (13): 

p = M[\^j)(^\} . (15) 

Substituting the Eq. (13) and taking the stochastic mean over J of the Gaus- 
sian correlation (14) we are led directly to the form (4) of IMS density op- 
erators. The 1-particle density matrix p and the parameter u are related to 
the correlation matrix a of the current by easily invertible matrix relations: 

a U P /-,^\ 
up = — — , a = . 16 

1 + a 1 — up 
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Measuring the 1-particle density matrix we could, up to the validity of the 
model, calculate the structure of the external current. Although such mea- 
surement is (so far) not completely possible we shall see in Sec. 5 that the 
inclusive correlation function gives the modulus of a directly. It is also seen 
from Eqs. (16) that a Gaussian shape, like 

Pv* ~ exp (-^k 2 + - l -R^_) , k + = ^, k_ = k-k', (17) 

for the 1-particle density matrix is not compatible with a Gaussian shaped 
current correlation matrix ctk'.k an d vice versa. 



5 Final state distribution vs. external current 



The final state distributions in IMS can be directly related to the currents 
J. The generating functional (5) can conveniently be re-expressed as an 
averaged functional over the fluctuating external current J: 



G[u] = M exp(^(u k -l)|J k | 2 j , 



(18) 



which is of course equivalent to Eqs. (6) or (9). 

The above equation has numerous useful consequences. The multiplicity 
distribution can be written in this form: 



p r = M[(]T|J| 2 ) r exp(-£|J| 



(19) 



while the factorial moments take the same form but without the exponential 
factor exp(—J2 \ J\ 2 )i i- e - : 



M[(£|J| 



(20) 



This phenomenon also characterizes the differences between the expressions 
of the exclusive and the inclusive distribution functions, respectively: 



/(I, 2, 



M 



I 7 |2 17-12 IT 



exp(-£|J|< 



(exclusive) 
(inclusive) 



(21) 
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as well as of the correlation functions. In particular, the inclusive correlation 
functions take the following form: 



C(l,2, 



M 



21 T \2 
Jl\ 



(22) 



the inclusive ones would contain the ominous exponential factor, too. The 
notation M[. . ] c means that in the 'expectation value' only the 'connected 
grafs' are to be taken into the account. In case of Eq. (22) it yields (r — 1)! 
'cycles', i.e. the 'cycle' ai2«23 ••• ovi and its variants for permutations of 
2,.,r| 

One can easily summarize the main result of this Section as follows. The 
counts nk, measured simultaneously in a collision event, are statistically iden- 
tical to the corresponding spectral intensities |Jk| 2 - Like their distributions, 
their corresponding moments are identical as well: 



(niTi2 • • • n r ) — M |Ji| 2 |J; 



\Jr\ 



(23) 



6 Bose-Einstein-condensation 

The IMS class of density operators (1) has a particular asymptotics. The 
'inverse temperature' matrix (3 must be positive. If it were degenerate the 
state (1) would not exist at all. A degenerate (3 can formally be interpreted 
as if the mode of zero eigenvalue became infinite hot. This mode is, in fact, 
becoming more and more populated but the infinite population is unattain- 
able. Nonetheless, an IMS with almost degenerate (3 would really be a Bose- 
Einstein condensate since this only requires a big finite number of bosons in 
a single quantum state. Speculations that the point of degeneracy, i.e. the 
point v = l/||p||, is the point of condensation (like in thermal Bose-systems) 
can not be verified for the IMS. 

Let us first recapitulate the kinematics of an IMS condensate. The con- 
densate mode does not interfere with the other modes so we can discuss it 
separately. We assume that our IMS is dominated by the condensate mode. 
The 1-particle density matrix has the form pk'k — WV 9 ! where <£>k is the 
condensate mode's wave function. If we introduce the condensate absorption 
operator a c = J^kfk^k then, using Eqs.(l-4), the condensate IMS can be 
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written as a thermal equilibrium state at temperature T = — 1/logi/: 

p c = (l-e-^)exp(-^ . (24) 

This state assumes a Hamiltonian a\a c which is not the real case, the con- 
densate is not even stationary in general. Yet, the form (24) is completely 
proper to calculate characteristics of the state by a thermal analogy. For 
instance, Eq. (5) yields directly the generating functional in the form: 

1 - e~ l ' T 

G[U] = l-exp(-l/T)E k ^ k | 2 ' (25) 
with the canonical thermal multiplicity distribution 

Pn = (1 - e-" T )e- n ' T (26) 

of mean multiplicity 

(»> = • ( 27 ) 

Let us observe that approaching the 'condensation point' corresponds to 
T — > oo and the population of the 'Bose-condensate' increases to the infinity 
while it is remaining thermally distributed all the time. 

Now I turn to the dynamic conditions for the fluctuating external current 
J. In the special case of the condensate IMS, the second relation in Eq. (16) 
becomes simply «k'k = (^)pk'k = (^Vk'Vk- Recall the definition (14) of a 
as the current's correlation matrix, which yields the following relation: 

M[J k ,J*] = (n)wtfi; • (28) 

Regarding that M[J^i should vanish by assumption (Sec. 4), the Gaussian 
fluctations satisfying the above relation must take the form 



■h = z^{n) Vk (29) 

for all k, where z is a random complex number of the standard Gauss- 
distribution (I/tt) exp(— \z\ 2 )d 2 z. Taking the stochastic mean of the modulus 
square of both sides we obtain: 

W = (n)\^\\ (30) 
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which also leads to 

("HEW- (si) 

k 

The Eqs. (29-31) show the simple way how the pulse of the effective current J 
determines the condensate wave function and the mean population. Actually, 
the mean multiplicity is identical to the overall intensity of the current pulse 
(31). The pulse's normalized spectral density is equal to the modulus square 
of the condensate wave function (30). The Eq. (29) seems, however, to be very 
restrictive since it imposes the same random phase and weight simultaneously 
for all current amplitudes Jk- 



7 Outlook: master equations 

I outline a possible generalization of the simple dynamics proposed in Sec. 4. 
Let us replace the Hamiltonian (12) by 

Hit) = g{t) (JZ(t)a k + J k (t)4) , (32) 

k 

where g(t) is a normalized function of characteristic width At, controlling the 
intensity of particle creation. The time-dependent currents Jk(i) are random 
functions of zero mean; let their correlator be of non-stationary white-noise 
type: 

g{t)M[Ju{t')Ji{t)\ = 6{t'-t)a k , k . (33) 

In the limit At —>■ 0, the random dynamics represented by Eqs. (32,33) 
reduces to the simplistic dynamics (12,14) of Sect. 4. The Hamiltonian (32) 
with the white-noise currents (33) yield the following master equation [OJ 
for the noise-averaged density operator (in interaction picture): 

-^ = 9^ "k'k (oifpait - ^{ak&k/, P}j ■ (34) 

One solves the master equation with vacuum initial state. If I /At is much 
greater than the typical pion energy then, via the relationships (16), the fi- 
nal state will tend to the IMS (1). For larger At the simple relations (16) 
do not hold. Though analytic calculations are still possible one can turn 
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to very powerful Monte-Carlo methods |jTB| developped for Markovian mas- 
ter equations. These MC algorithms will yield the density operators of the 
multiparticle final states without " struggling through" |TJ| the usual Wigner- 
function formalism. 



8 Summary 

The aim of the Letter was partly pedagogical. To avoid Bose-symmetrization 
'by hand', I used standard quantum mechanical considerations to construct 
and to analyze the 'independent multiboson states'. I showed how these 
states emerge from a simplistic version of chaotic current models and I derived 
the relationship between the IMS states and the correlator of the currents. I 
briefly recapitulated the generating functional representation of multiparti- 
cle counts. Beyond methodological matters, I found that the 'Bose-Einstein 
condensate' would be thermally populated and the condensation point cor- 
responds to the infinite hot state. I restricted my analysis for the IMS of 
Sec.l and for the simplest chaotic current mechanisms with trivial time- 
dependences. The simple choice allows for transparent relationships between 
the current and the final multiparticle state. (Other works, like e.g. Ref . [|1C|1 . 
incorporate a more realistic space-time evolution of the multiparticle source.) 
There is, nonetheless, a particular advantage of any underlying dynamics 
whether realistic or not. Usually it allows economic simulation methods for 
the physical quantities of interest. To this end, I proposed a quantum master 
equation known to be suitable to efficient Monte-Carlo simulations. 
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